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Abstract 

The mean density of a random closed set Q in R d with Hausdorff dimension n is the Radon- 
Nikodym derivative of the expected measure E[7i n (S n ■ )] induced by O with respect to the 
usual d- dimensional Lebesgue measure. We consider here inhomogeneous Boolean models 
with lower dimensional typical grain. Under general regularity assumptions on the typical 
grain, related to the existence of its Minkowski content, and on the intensity measure of the 
underlying Poisson point process, we prove an explicit formula for the mean density. The 
proof of such formula provides as by-product estimators for the mean density in terms of the 
empirical capacity functional, which turns to be closely related to the well known random 
variable density estimation by histograms in the extreme case n — 0. Particular cases and 
examples are also discussed. 
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1 Introduction 

In many real applications it is of interest to study random closed sets at different Hausdorff di- 
mensions and their induced random measure; in particular several problems are related to the 
estimation of the density, said mean density, of the expected measure 

E[/x e ] (B) := E[H n (0 n B)} \/B 6 £ Rd 

induced by a random closed set : (f2, J, P) — > (F, cry) in M. d with Hausdorff dimension n 
(H n , B^d, F and <tf denote here the n-dimensional Hausdorff measure, the Borel cr-algebra of 
M. d , the class of the closed subsets in M. d and the cr-algebra generated by the so-called hit-or-miss 
topology [16] , respectively). While the extreme cases n — d and n = are easy to handle with 
elementary analytical tools, problems arise when < n < d. Indeed, while the expected mea- 
sure induced by a d-dimensional random set is always absolutely continuous with respect to the 
d-dimensional Lebesgue measure and its density can be easily obtained by applying Fubini's theo- 
rem (as stated in the early work by Robbins |19j). and the mean density of a random point is given 
by its probability density function (and so the problem of its estimation has been largely solved 
since long in nowadays standard literature, by means of either histograms, or kernel estimators 
(e.g. see |21))V when we deal with a general lower dimensional random closed set it can be more 
demanding to check that the induced expected measure is absolutely continuous and to compute 
and estimate its density. Even if inhomogeneous random closed sets appear frequently in real 
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applications, only the stationary case has been extensively studied so far (e.g., see [71(22] and ref- 
erence therein). Indeed, although some results about the mean densities of certain inhomogeneous 
random closed sets are available in current literature (e.g., see [5J [T3J [T^l [201 121] ) > mainly via tools 
from integral geometry and stereology, the pointwise estimation of the mean density of nontrivial 
lower dimensional inhomogeneous random closed sets is still an open problem, whenever gradient 
structures [13 or local stationarity cannot be assumed. 

Our goal is to obtain a pointwise result for the mean density of a wide class of inhomogeneous 
random sets which might be of interest from a statistical point of view as well. We consider here in- 
homogeneous Boolean models with lower dimensional typical grain, being considered basic random 
set models in stochastic geometry [3] , but we do not exclude that our approach could be applied to 
other kinds of inhomogeneous lower dimensional random closed sets in further developments. The 
main result of the paper is stated in Theorem where we prove an explicit formula for the mean 
density; we point out that, even if such formula could be obtained in a more direct way via the well 
known Campbell's formula (see Remark Hip , the proof we propose here provides as a by-product 
estimators of the mean density, answering, in the case of Boolean models, to an open problem in 
PP. To this end, the classical Minkowski content notion plays a central role throughout the paper; 
as a matter of fact, the regularity assumptions on the typical grain we require are closely related 
to the well known general assumptions which guarantee the existence of the Minkowski content of 
deterministic closed sets. 

The paper is organized as follows. In Section[2]we recall some basic definitions and preliminary 
results useful for the sequel. In Section [3] we prove a generalization of the Minkowski content 
of closed subsets of M. d (Theorem [6]); we will apply such result in the final part of the proof of 
Theorem [9l In Section [4] we introduce a wide class of Boolean models in K d with lower dimensional 
typical grain. We prove that the expected measures induced by the Boolean models we consider 
are absolutely continuous with respect to the usual Lebesgue measure fLemmallOp and we provide 
an explicit formula for their density. Such formula simplifies in the special cases in which the 
Boolean model is assumed to be stationary ( Corollary I12[) . or to have deterministic typical grain 
fCorollary|13p. A simple example of inhomogeneous segment Boolean model and links with current 
literature are also provided. Finally, the problem of the estimation of the mean density is considered 
in Section [51 in particular, we define an estimator of the mean density (Proposition I15p . which can 
be considered as the generalization to the case < n < d of the classical random variable density 
estimation by histograms in the extreme case n = 0. 

2 Basic notation and preliminaries 

We recall that TC° is the usual counting measure, H n (B) coincides with the classical n-dimensional 
measure of B for 1 < n < d integer if B € Z? R d is contained in a C 1 n-dimensional manifold 
embedded in M. d , TL d (B) coincides with the usual ei-dimcnsional Lebesgue measure of B for any 
Borel set B C R d . Throughout the paper dx stands for Ti d (dx). A closed subset S of M. d is said to 
be countably Tt n -rectifiable if there exist countably many n-dimensional Lipschitz graphs Yi C K d 
such that S \ UiTi is 7Y ra -negligible. (For definitions and basic properties of Hausdorff measure and 
rectifiable sets see, e.g., [31 [TUl [H] -) We call Radon measure in R d any nonnegative and cr-additive 
set function /i defined on B R d which is finite on bounded sets, and we write /x <C TL n to say that 
/i is absolutely continuous with respect to TL n . We will say that a random closed set O satisfies a 
certain property (e.g. has Hausdorff dimension n) if Q(lu) satisfies that property for P-a.e. ui £ Q. 
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Throughout the paper we will consider countably 7i™-rectifiable random closed sets in M. d , with 
1 < n < d — 1 integer, such that E[/ie] is a Radon measure; the particular cases n — and n = d 
are trivial. (For a discussion about measurability of Tt n (Q) we refer to [2"5I [5].) 
Whenever E[/^e] is absolutely continuous with respect to H d , the following definition is given [51125] 

Definition 1 (Absolute continuity in mean and mean density) Let be a countably H n - 
rectifiable random closed set in R d such that K[fj,&] is a Radon measure. We say that is absolutely 
continuous in mean z/E[^e] <C TC d . In this case we call mean density of Q, and denote by Xq, the 
Radon-Nikodym derivative o/E[/ie] with respect to 7i d . 

Remark 2 fn the case n = with = X random point in R d , we have that is absolutely 
continuous in mean if and only if X admits a probability density function fx, and so \® = fx- 
On the other hand, it is easy to see (as an application of Fubini's theorem in f! x R d with the 
product measure ¥xTl d ) that any e?-dimensional random closed set in R d is absolutely continuous 
in mean with mean density \q(x) = P(x £ 0) for H d -&.e x £ R d . 

The problem of the approximation of the mean densities in the general setting of spatially inho- 
mogeneous processes has been recently faced in pQ, where an approximation, in weak form, of the 
mean density for sufficiently regular random closed sets is given in terms of their d-dimensional en- 
largement by Minkowski addition. More precisely, denoting by 5® r := {x £ M. d : 3y £ S with \x — y\ < r} 
the closed r-neighborhood of a closed set S C R d , and by B r (x) the closed ball centered in x with 
radius r, for any compact window W C R d let rV(0) : — > M be the function so defined: 

Tiy(0) := sup{7 > : 3 a probability measure r\ <C H n such that 

r](B r (x)) > 7 r w VxeenW ffil , r £ (0,1)}; (1) 

then the following theorem holds [T]. 

Theorem 3 Let be a countably TC n -rectifiable random closed set in M. d such that E[/xe] is a 
Radon measure. Assume that for any compact window W C M. d there exists a random variable Y 
with E[y] < oo, such that l/r^(0) < Y almost surely. If is absolutely continuous in mean, 
then 

lim / — j S "' dx — / Ae(^)dx (2) 

rloJ A b d _ n r d - n J a 

for any bounded Borel set A C M. d with H d (dA) = 0. 

Hence the above theorem gives a weak result for the mean density of very general lower dimensional 
random closed sets. In order to obtain a pointwise result, using the fact that A is arbitrary, we 
should prove that limit and integral in @ can be exchanged; in such a way we could state that 

A e f x ) = Hm € ^ 0r) for H d -&.c. x £ R d . 

V ' riO b d ^ n r d - n 

and, as by-product, the right side of the above equation could suggest estimators for Xq(x) in 
terms of the capacity functional of 0. The proof of the validity of this formula for absolutely 
continuous in mean random sets seems to be a quite delicate problem, with the only exception 
of the stationary ones and the extreme cases n — d and n = 0. In Theorem [S] we prove that ([3]) 
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holds for a wide class of (inhomogeneous) Boolean models and, in particular, we give an explicit 
representation of the mean density Ae in terms of the intensity and of the typical grain of the 
process. To this end we recall now the definition (typically standard) of Boolean model and we 
refer to [H [S] for details about point processes and related concepts and results. 

Definition 4 (Boolean model) Let = {xijigN be Poisson point process in R d with intensity 
f and let {Zj}j g pj be a sequence of i.i.d. random compact sets in M. d , which are also independent 
of the Poisson process $ . Denoting by Zq a further random compact set of the same distribution 
as the Zi 's and independent of both of them and of the resulting random set 

9 := \J( Xi + 

i 

is said (inhomogeneous) Boolean model with intensity / and typical grain Zq. 

In our assumptions Zq will be a lower dimensional random closed set uniquely determined by a 
random quantity in a suitable mark space K, so that Zq{s) is a compact subset of M. d containing 
the origin for any seK. We remind that in common literature is usually assumed that 

E[card{i : (x. t + Z t ) n K ^ 0}] < oo V compact K C R d (4) 

(card stands for cardinality), and that a Boolean model as above can be described by a Poisson 
point process in R d x K with intensity measure A(d(cc, s)) = f(x)dxQ(ds), where Q is a probability 
measure on K representing the distribution of the typical grain. 



3 A generalization of the Minkowski content of closed sets 

Denoted by bk the volume of the unit ball in R fc , the n-dimensional Minkowski content A4 n (S) of 
a closed set S C R d is defined by 

M n (S) := lim H<i{S T } 

rlQ b d _ n r d - n 

whenever the limit exists finite. 

General results about the existence of the Minkowski content of closed subsets in R d are known 
in literature, related to rectifiability properties of the involved sets. In particular, the following 
theorem is proved in [3j (p. 110). 

Theorem 5 Let S C M. d be a countably 7i n -rectifiable compact set and assume that 

r](B r (x)) > 7 r™ Va; € S, Vr e (0, 1) (5) 

holds for some 7 > and some Radon measure r\ in R d absolutely continuous with respect to 7i n . 
ThenM n {S) = H n {S). 

Note that such theorem extends the well-known Federer's result ([11] . p. 275) about the existence 
of Ai n (S) for n-rectifiable compact sets S C R d (i.e. S is representable as the image of a compact 
subset of K ra by a Lipschitz function from R" to R d ) to countably 7i"-rectifiable compact sets (see 
Remark 2.3 in [5]). Moreover, in many applications condition §5$ is satisfied with r](-) = H n (S n • ) 
for some closed set S D S, and for n = d— 1, it is not hard to check that such condition is satisfied 
by all sets with Lipschitz boundary (see [HI3])- 
We state now the main result of this section. 
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Theorem 6 Let fx be a positive measure in M. d absolutely continuous with respect to 7i d with 
density f such that 

i) f is locally bounded (i.e. snp x£K f(x) < oo for any compact K C M ); 

ii) the set of all discontinuity points of f is Ti n -negligible. 

Let S C M. d be a countably 7i n -rectifiable compact set such that condition ([5]) holds for some 7 > 
and some probability measure r\ in K d absolutely continuous with respect to Tt n . Then 

rio b d - n r d n J s 

Remark 7 The above theorem may be seen as a generalization of TheoremjS] indeed, the classical 
Minkowski content follows as particular case by choosing / = 1 and noticing that if a Radon 
measure 77 as in Theorem [5] exists, then it can be assumed to be a probability measure without loss 
of generality. Indeed, it is sufficient to consider the measure rj(-) := n(W n ■ )/rj(W), where W is a 
compact subset of M. d such that S®i C W; it is clear that fj is a probability measure satisfying 

v(B r (x)) > VxeS, Vr e(0,l). 

Furthermore, a classical result from geometric measure theory (e.g., see [3] Theorem 2.56) tells 
us that if [i is a positive Radon measure on R d and B 6 B^d such that 

limsu P M( f r(a:)) >te (0,oo) VxeB, 
rio Kr n 

then /j(-) > tTL n {B n •). Hence, any set 5* C K d as in Theorem[H]has finite 7Y"-measure. 

In order to make the proof of Theorem [S] more readable, we remind that Lemma 15 and Lemma 
16 in [lj provide a local version of Theorem [5] and an upper bound for the Minkowski content of 
compact sets in R d , respectively; for our purpose we summarize as follows. 

Lemma 8 If S C M. d is a countably TL n -rectifiable compact set such that condition ((SJ) holds for 
some 7 > and some finite measure 77 in M. d absolutely continuous with respect to 7i n , then 

KHM < W < 2, (6) 

b d - n r d - n ~ 7 b d -n 

and 

ny^nA) =H n {SnA) (7) 

rio b d - n r d - n K ' K ' 

for any A e B R a such that H n (S D dA) = 0. 

Proof of Theorem® It is well known that, since / > 0, there exists an increasing sequence {fk}ken 
of step functions 

N(k) 

f k ( x ) - J2 4 fe)l Af » (*)> a ? ^ °' Ar ( fc ) e N > 

3=1 
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converging to /. 

By ii), such sequence {/&} can be chosen so that 7i n (S n dAj) = for all j, k. 
Let 

*< r > = -fc*^ 

We may observe that 

(a) liWr) = £ «f to ".I? } § E ^^(5 n Af ) = / f k ( X )H"(dx). 

riO f J riO Od- n r a " "T-f Js 

Since S is compact with H n (S) < oo by Remark [3 and / is locally bounded, it follows that 
/ fe (x)H"(dx) < oo. 

J s f(x)dx 

(b) lim gk(r) = gr — -5 uniformly in (0, 1); indeed: 

fc^oo b d ^ n r d n 

we have that /fc | / uniformly in S*0. r for any r > because / is bounded in 5*© r , being S® r 
compact, i.e. for all e > there exists kg such that 

sup \f k (x) -f(x)\ < e Vk>k . (8) 
Hence, for any e > 0, for all k > ko 



9k(r) ~ -7 -J-- 



< J s JMx)- f(x)\dx o o i ^ ^ 

b d - n r d - n b d _ n r d - n ~ 7 6 d _„ 

As a consequence of (a) and (b) lim r jo hm^—too gk{r) — limfc^oo lim r j.o gk( r )- Thus the following 
chain of equalities holds 

,. n(S@r) ,. /s fflP /( a, ) dB ,. / Sffir limfc^oo /fc(x) da; fs^M*) 6 * 
lim j — = hm — r 2 j = lim — ; r = hm hm — ^ 3 

''10 b d _ n r d - n rio b d - n r d - n no b d _ n r d - n r|0 fc^oc 6 d _„r d " n 

Js fflr (a) 



= lim lim Js ®- ** ( ^ lim / / fc (z)H"(dx) = / /W(dx). 



□ 



4 Mean density of inhomogeneous Boolean models 

Let the Boolean model so defined 

0(w) := U + Z (ai) e (I, 

(ai,Si)S#(w) 

with $ Poisson process in R d x K and Zo typical grain in R d , satisfying the usual condition Q 
for Boolean models. To lighten the notations, from now on we denote by 

- Z x > s :=x- Z {s) V(i, s) G R d x K; 

- i5 the random variable on K so defined 5(s) := diamZo(s), where diam stands for diameter; 

- Eq the expectation with respect to the probability measure QonK. 



G 



Note that the condition (0J is equivalent to 

E[$({(y, s) £ R d x K : (y + Z Q (s)) n B R (0) + 0})] < oo Vi? > 0, 
and so, in terms of the intensity measure A of $, to 

A(dy x ds) < oo Vi? > 0, (9) 



where — Zq is the symmetric of Zq with respect to the origin. 

Assumptions: Let $ have intensity measure A(dy x ds) = f(y)dyQ(ds) satisfying (J9j> ; further, 
let us assume that the following conditions on / and Zq are fulfilled: 

(Al) Zq(s) is a countably 7i"-rectifiable compact set in R d for Q-a.e. s 6 K. Further there exist 
7 > and a random closed set Zq D Zq with Eq[H u (Zq)) < oo such that, for Q-a.e. s € K, 

H n (Z (s) n B r (x)) > 7 r" Vx G Zq(s), Vr e (0, 1). (10) 

(A2) the set of all discontinuity points of / is ^"-negligible and / is locally bounded such that 
for any compact set K C M. d 

sup f{y)<U (11) 
holds for some random variable £k with Eq[H'' 1 (Zq)£ : k} < oo. 

Theorem 9 (Main result) Any Boolean model as in Assumptions is absolutely continuous in 
mean with mean density Xq given by 

A e (x)= / f f(y)H n (dy)Q(ds) forH d -a.e.xeR d . (12) 

Before proving the above theorem, let us make a few observations about the above Assumptions 
in order to clarify their level of generality. 

• Condition is trivially satisfied whenever / is bounded, or / is locally bounded and 
diamZo <c£ K + Q- almost surely. 

• Assumption (Al) is satisfied by a great deal of typical grains Zq with Zq = Zq or Zq = ZqUA, 
with A sufficiently regular random closed set to control the case when Zq can be arbitrarily 
small (see [3 [23]), and it assures that Zq is a countably 7i™-rectifiable random compact set 
with finite expected 7i"-measure. 

• Stationary case. If O is stationary with / = c > 0, then only the regularity assumption 
(Al) on the typical grain Zq is required. Indeed, the assumption (A2) is trivial and the usual 
condition ©, which in this case becomes 

E Q [H d ((Z ) eR )} <oo Vi?>0, 

is satisfied thanks to the assumption (Al): for all R < 2 by Lemma|S]with n(-) = Tt n (ZQ n • ) 
we have that E Q [n d {{Zo) mR )] < E[H n (Z )]2 n 4 d b d R d - n / 7 , while for all R > 2, by repeat- 
ing the same argument of the proof of the quoted proposition, it is easy to check that 
E Q [H d ((Z ) BR )] < E[n n (Zo)}2 n (AR) d b d /-f. 
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• Deterministic typical grain. If has deterministic typical grain Z$ C K d the above 
Assumptions can be replaced by: has intensity / and typical grain Z§ = S as in the 
hypotheses of Theorem [5] (See Corollary Q]|] for details.) 

Lemma 10 For any Boolean model as in Assumptions, E[/ie] * s a locally finite measure abso- 
lutely continuous with respect to 7i d . 

Proof. For any R > let B R := {(y, s) e M d x K : (y + Z (s)) D B R (0) ^ 0}; then 

E[H n (Q n B R (Q))] = E[E[?T(0 D Bfl(O)) | < E Q [^"(Z )]E[$(S i? )] < oo 

by Q and condition (Al); so E[/te] is locally finite. 

By contradiction, let E[/ie] be not absolutely continuous with respect to H d ; then there exists 
A C M d with H d (A) = such that E[W n (0 n A)] > 0. In particular, 

o<p(H"(enA)>o)<p( w n ((xi + Zo{si))r\A) >o) = p(«(.4) >o), 

(xi,si)e* 

where 

.4 := {(y, s) £ K d x K : + Z (s)) nA)> 0)}. 

Denoting by „4 S := {y G M d : (y, s) e .4} the section of i at s e K, we may apply Fubini's 
theorem to get 

fn n ((y+Z Q (s))nA)dy=f(f l A _ y (x)H n (dx))dy = f (f l A _ x {y)dy)H n (dx) = 0, 

JA, ./.A s K Jz (s) ' JZ a (s) K JA S ' 

because H d (A) = 0. Being the function y i— > H n ((y + Z (s)) n A) strictly positive in „4 S , we 
conclude that H d (4 s ) = for all s <E K. Then it follows 

E[9(A)] = J A A (dy x da) = J (J /(t/)dy)Q(d S ) = 0; 

but this is impossible, because P($(.A) > 0) > implies E[$(.4)] > 0. □ 

Proof of Theorem fp) Clearly is a countably 7Y"-rectifiable random closed set in R d , and by 
Lemma [TU] it is absolutely continuous in mean, so E[/ie] = Ae7i d for some integrable function Ae 
on R d . 

Let W be a fixed compact subset of R d . For any w s 51 let us consider the set 

0(w) := (J Xi + Z (si) 

and the probability measure r/yy absolutely continuous with respect to TL n so defined 

r(enn%nfl) 

vwyB) ■= ~ vBeB R d. 

w»(e(w)nw ffl i) 

Note that for any a; £ 0(w) (1 W®i there exists (x, s) £ $(w) such that i£x + ^o(s), and so 
MiW) > 9 — ^ r « W € (0, 1). 

H"(e( w )n%) " w(e(w) n w e i) 
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Thus, here the function rV(0) defined in |T]) is such that r^(6) > 7 /TL n (Q (~l by the same 

argument in the proof of Lemma ITD1 it is easy to check that E[H"(8nWgi)] < 00, then Theorem[3] 
applies and we get 

E ^ ) = S/, ! B^ d * (13) 

for any bounded set A G B^d such that H d (dA) = 0. 
Let us denote by Z x,r the subset of R d x K so defined 

Z x > r := {(y, s ) G R d x K : x G (y + Z (s)) e r} - {(y, s) G R d x K : y G Zg}, 

and observe that Vi G A, Vr < 2 

p(x g e ffi r) _ P($(^ x ' r ) > 0) 1 - e - A ( 2X,r ) 



b d - n r d - n bd- n r d - n b d - n r d - 

A(Z x - r ) _ 1 
bd- n r d ~ n ~ b d - n r 



(14) 



< r^^ = 7— ^r/ / /(v)dvQ(d*) 



< / -p /(V)0(d*) 

GD,GD 2 n A d bd f ~ (-A2) 

< / W"(^ (s))^ ffi2 ( S Q(d S ) < oo. 

Then by the dominated convergence theorem we can exchange limit and integral in Eq. (fT5| . 
Similarly, we have that for all r < 2 

! f f[ f v < «y*Mw> »p /<») < ^«"(2o (8 ))f(») 

r d-n — U r d-n Je a ~ rvhj 

Ud—n 1 Ud—n 1 v pzZ. I u d—n 

for some random variable £ with Eq[7Y™(Zo)£] < 00, by © and (A2). Then, the dominated 
convergence theorem implies that 

lim , 1 J n = / lim -Si d _ K Q(ds) = / / f(y)H n (dy) Q(ds), (15) 

where the last equality follows by applying Theorem [H] with /i = fH d and S = Z x,s . 
Summarizing, we have that for any bounded set A G B R d with Tt d (dA) = 

E u e A = hm / -H ^dx= / hm — ^-dx 

rio J A b d - n r d - n J A rio b d - n r d ~ n 

m f ,. l-e- A ( z *' r ) , /• A(Z K ' r ) + o(r d -™) , 
= / hm — ; r dx = hm ; dx 

7^ r;o b d - n r d - n J A HO b d - n r d - n 



/(y)W"(dy)Q(d S ) dx 



KJZ 



We conclude that E[/te] has density 



Ae(a;) 



/ / f {y)H n {dy)Q(ds) for H d -a.e. x G 



□ 
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Remark 11 Formula (fT2")) could also be obtained in a more direct way by the well-known Camp- 
bell's formula (e.g., see [3J, p. 28, and [T3] for a similar application), after having shown that 
E[/ie](-) = E[^ i 7<"((xj + Zj) n • )] for any Boolean model 9 as in Assumptions, being zero the 
probability that different grains overlap in a subset of M. d of positive 7i n -measure. (To show such 
property, one could proceed arguing similarly to the final part of the proof of Lemma 1101 ) On 
the other hand, a proof via Campbell's formula seems not to give any hint to provide computable 
estimators for the mean density (see Section [5] for a more detailed discussion) . 

In the particular cases in which is stationary, or the typical grain Zq is a deterministic subset 
of M. d satisfying the hypotheses of Thcorcm[5](with r\ probability measure by Remark^, Thcorem[9] 
specializes as follows. 

Corollary 12 (Stationary case) Let be a stationary Boolean model with intensity f = c > 
and typical grain Zq satisfying the assumption (Al). Then is absolutely continuous in mean with 
mean density 

Ae(a:) = cE Q [H n (Z )} Vx e K d . 

Proof. We have observed that satisfies the Assumptions and Ae is constant on R d . Then the 
thesis follows by TheoremEE noticing that H n {Z x <") = H n (Z (s)) for all (a;, s) el^xK and that 
Ae is constant being E[/xe] translation invariant. □ 



Corollary 13 (Deterministic typical grain) Let be a Boolean model in R d with determinis- 
tic typical grain Zq C M. d . If the underlying Poisson process $ in M. d has intensity f locally bounded 
and such that the set of all its discontinuity points is TL n -negligible, and if Zq is a countably 7i n - 
rectifiable compact set such that condition ^ holds for some 7 > and some probability measure 
7] in M. d absolutely continuous with respect to 7i n , then is absolutely continuous in mean with 
mean density 



Xe(x)= f(x-y)H n (dy) for H d -a.e. x el 4 . 

Jz 

Proof. Since Zq is compact diamZo = S < 00, and we know by Remark [7] that TL u {Zq) is finite. 
Then the assumption (A2) and the condition §§§ are easily checked. It follows that E[/ie] is locally 
finite and, by proceeding as in Lemma [TOl again we have that E[/ze] "C 'H d . 
For any y £ K d let rj v be the measure on M. d so defined 

rf{B):=n{B-y) MB £ B Rd . 

Let W be a fixed compact subset of R d and for any ui £ consider the measure 

V{B > '- card{x. ( i $M : (x. t + Zq) H W m + 0} V ^ e 

Note that 

- 77 is a probability measure absolutely continuous with respect to 7i n ; 

- for any x £ 0(o>) PI W^i there exists x £ $(w) such that x £ x + Z , and so 

^ Br{x)) ~ cardie $H : (J, + Zq) n W el P ®f * G ( °' 1} ' 
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Then 9 satisfies the hypotheses of Theorem[3] with F = cardjxj G $(w) : (xi + Zo)nW^i 7^ 0}/7, 
which has finite expectation thanks to ©, and so again we have that Eq. (|13[) holds for any 
bounded set A G B R d such that Ti d (dA) — 0. Finally, by proceeding as in the final part of the 
proof of Theorem [5] with deterministic Z , we conclude that still holds and 

A (x) = / f(y)H n (dy) = f f(x- y)H n (dy) for H d -z.e. x G R d . 

J'X-Zq J Zq 

□ 

Remark 14 (The special case n = d — 1) A well known tool in stochastic geometry is the so- 
called local spherical contact distribution function, defined as i?e( r j %) '■— 6 @©r \ x £ @) for all 
(r, x) e R+ x R d (e.g, see [H]). If 9 is a Boolean model as in Assumptions with (d— l)-dimensional 
typical grain, then P(x G 9) = for H d -&.e. x G R d , and so the mean density Ae can be written 
in terms of Hq as well 

M*0 - l-w-H e (r,x) lr=0 for W d -a.e. x G R d . 
2 or 1 

Therefore known results about Hq could be applied to (d — l)-dimcnsional Boolean models to 
obtain further information on their mean density. 

We conclude this section by giving a very simple example of inhomogeneous Boolean model in 
R d , in order to show how the mean density can be easily computed by applying Theorem 

Example 1 (Segment Boolean model) For sake of simplicity we consider a Boolean model 9 
of segments in R 2 , but a similar example can be done in R d with d > 2. 

So, let K = R + x [0, 2tt] and for all s = (I, a) G K let Zq(s) be the segment with length I and 
orientation a so defined 

Zq(s) :— {(u, v) G R 2 : u = t cos a, v = t sin a, r G [0, 1]}. 

We consider the case in which both length and orientation are random. Denoting by L the Re- 
valued random variable representing the length of Zq and by Pi(dZ) its probability law, let $ be 
the marked Poisson process in R d x K having intensity measure A(dy x ds) = f(y)dyQ(ds) with 
v) = u 2 + v 2 and Q(ds) = -^daF^dl) such that J R+ l 3 ¥i,(dl) < 00. (This last assumption is 
to guarantee that the usual condition ^ is satisfied; for a different intensity / we might have a 
different condition on the moments of L.) It is easily shown [UH3] that 9 satisfies the assumption 
(Al); by noticing that 

I f{y)dy<{l + R) 2 {2lR + TTR 2 ) Vs = (I,a)eK, 

J(-Z ( S )) @R 

it follows that 

f ( A(dyxds)< f (l + R) 2 (2lR + TTR 2 )P L (dl) < 00 Vi? > 0, 

JkJ(-Z„(s)) 9r Jr + 

so condition |9]) (and similarly the assumption (A2)) is satisfied. Hence Theorem |9] applies and we 



11 



got 



\e(xi,x 2 ) = 



poo 2 p2tt rl 

/ — / / f(xi — rcosa, X2 — t sina)dr daPi,(dZ) 
Jo 2vr J J 

/■°° i /" 27r / i \ 

/ — / ( ( x i + x l) 1 ~ ( x i cosa + x 2 sina)/ 2 + -I 3 )daP L (dl) 
Jo 27r J V 3 / 



= {xj + xl)E[L] + -E[L 3 }. 

3 

Note that in the particular case / = c > 0, is stationary and by Corollary [12] we obtain the 
well known result X e (x) = cE[L] \/x G R d (cf. [7j,p.42). 



5 On the estimation of the mean density 

While the problem of the mean density estimation has been widely examined in the stationary case, 
the inhomogeneous one has been mainly faced by assuming local stationarity. For instance, is 
assumed to have a gradient structure, i.e. it is considered to be homogeneous perpendicularly to a 
particular gradient direction (see |13j). or it is assumed to be homogeneous in certain subregions 
of R d , so that the known results in the homogeneous case can be applied to estimate a stepwise 
approximation of the mean density Ae- We also mention that the stationary Boolean model is 
often assumed to have unknown constant intensity c > and known mark distribution Q, so that 
only the parameter c has to be estimated, being Ae = cEq[T-C 1 (Zo)] in this case. A series of results 
about the estimation of the intensity c of the underlying Poisson point process associated to 0, 
related to the estimation of Ae, can be found in [7] §3.4 (see also [HI [22]). 

Having now an explicit formula for the mean density of inhomogeneous Boolean models as in 
Assumptions, in all situations in which it is possible to estimate the intensity / and the mark 
distribution Q of the typical grain Zo, an estimation of Ae could be obtained. Actually, the 
estimation of / and Q as well as the computing of the mean density Ae (x) at a given point x £ M. d 
might be quite hard. In this section we introduce an estimator for the mean density Ae(a^) of 
Boolean models as in Assumptions in the general case of / and Q unknown. 

In the proof of Theorem [5] we have shown, in particular, that Ae is also given by the limit in 
© for H d -&.e. x € R d . Noticing that P(x £ ffi r) = T @ (B r (x)), where T e is the capacity (or 
hitting) functional of [16] , a natural estimator of Ae (x) can be given in terms of the empirical 
capacity functional of 0, without estimating the intensity / and the distribution of Zo separately. 
We recall that the empirical capacity functional T~ based on an i.i.d. random sample 3j,, . . . , Sjy 
of a random closed set 2 is defined as (see, e.g., jT^] ) 

1 N 

T£(K) := - ls,n^0, V compact K C M d , 
i=i 

and that the strong law of large numbers implies that Tff(K) converges almost surely to T~,(K) 
for any compact subset K of W l . 

Let be an inhomogeneous Boolean model in M. d as in Assumptions and 0i, . . . , 0jy be an i.i.d. 
random sample of 0; for any fixed r > we have that 



Tg(B r (x)) ] _ P(x e e r) / TgjBrjx)) ^ _ P(x£9 $r )(l-P(x£8 $r )) 



bd- n r 



d—n 



b d - n r d - n ' V b d - n r d - n J N{b d - n r d - n f 1 ' 
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Hence, © and p^|) suggest that we take 



with itV such that 



lim i? w = and lim NRt n = oo, (18) 

N~>oo N->oc " 

as asymptotically unbiased and consistent estimator of the mean density Ae(x) of at point x. 



Proposition 15 Let O be a Boolean model in M. d as in Assumptions and {Oi}igN be a sequence 
of random closed sets i.i.d. as 0; then 

lim Aq(x) = \q(x) in probability, for Tt d -a.e. x G 

A r ^oo 

Proof. By the definition of Aq it follows 



^ = Aq(x) tor 7t -a.e. a; G 



lim E[AS(a:)]= Urn v J^ wy = A e (x) 

We check now that the variance of Aq (x) goes to 0. Since 

lim P(x G e efljv ) =P(ie6)=0 for H d -a.e. because n < d, 

iV— »oc 



>(x € e ffifl J 



J im "I Jd-n = X ®( x ) e M for ^ d -a-e. a; £ R d by ©, 



and 



AT 



A^oo NR 



we have that for H d -&.e. x G M d 

urn var Ao u = hm -, = 0. 

Hence the thesis follows. □ 

Then, a problem of statistical interest could be to find the optimal width Rn satisfying condition 
(fT5|) which minimizes the mean squared error of Aq (x) (i.e. E[(A@(x) — Aq(x)) 2 ]). To investigate 
this problem is not the aim of the present paper and we leave this as open problem for further 
developments; we point out here that Aq can be seen as the generalization to the case of n- 
dimensional random closed sets, of the well known estimator of the probability density of a random 
point, which is a particular 0-dimensional random closed set. 

Remark 16 (The special case n — 0) Even if the particular case n = can be handle with 
much more elementary tools, it is easy to check that if = X is a random point in M. d with 
probability density function fx, Eq. §S§ holds with Xx = fx (it is sufficient to observe that 
E[H°(X n ■ )] = ¥(X G • )) and the estimator A^ turns to be closely related to the well known 
definition of histogram (see [23 for details). Let us consider the case of a random variable X with 
density fx', then Proposition 1151 applies with d = 1 and n — 0, making explicit the correspondence 
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with the usual density estimation by means of histograms (e.g., see |18j §VII.13). Indeed, if 
is a sequence of i.i.d. random variables with the same distribution of X , we define 

t , x 'Tn i s GD Sill ^Br {x)[Xi) card{i : X, £ I x } 
f X (x) := \ x (x) = = ^ , 

where I x is the interval in ]R centered in x with length \I X \ = 2i?jv with the usual condition 

\I X \ — > and N\I X \ — > oo as N — > oo. 

Therefore, statistical problems and techniques related to the choice of the "optimal width" Rjy in 
the general case N > could be investigated starting from available results for random variables. 

We conclude observing how the following corollary of Theorem [§1 which could lead to further 
developments on this topic and on the estimation of the mean density, is consistent with a couple 
of available results in literature. 

Corollary 17 Any Boolean model as in Assumptions is absolutely continuous in mean with 
mean density 

* M r IE[card{(x. t , 3i ) £ $ : + Z (s t )) n B r (x) £ 0}] , , 

\q(x) — lim -, for n -a.e. x £ K . (19) 

rj.0 b d - n r d ~ n 

Proof. The assertion follows directly by (|T5|) . noticing that A(Z x - r ) = E[ca,rd{(xi, si) G $ : 
(s i + Zo(« i ))nfl r (aj)^0}]. □ 

By Proposition 21 in [I] we get that for a locally finite union O of i.i.d. random closed sets Ei 
with Hausdorff dimension n < d it holds 

P(s g e ffir ) = EtcardjE ^n^)^}] for ^ Rd 
rio b d - n r d ' n no 6rf-„r d -" v ; 

provided that at least one of the two limits exists. On the other hand, in [2U] the mean density of a 
class of nonstationary n-fiat processes is studied and a similar result to (|19p is obtained. Namely, 
we remind that a n-flat process in M d (with n integer less than d) is a point process X on £„, 
the space of n-dimensional planes in M. d ; it is proved that if the intensity measure of X has a 
continuous density h with respect to some translation- invariant, locally finite measure on £ d , then 
the n-dimensional random closed set G := Usex ^ nas continuous mean density 



\ e (x) = / h(x + £)*(di), 

where £^ is the Grassmannian of n-dimensional linear subspaces in M. d and \& is a finite measure 
on coming from a decomposition result of the intensity measure of X, and in particular it is 
claimed that (see [H], p. 142, or [4 , p. 179) 

, , , E[card{£: e X : E n BJx) + 0}] 
Ae (x J = hm - . 

Ho b d _ n r d ~ n 

We may like to notice that Theorem [3] applies to n-flat process; hence we are lead to conjecture 
that the exchange between limit and integral in ^ may hold for further processes = [j-Ef, 
union of i.i.d. n-dimensional random closed sets, so that Proposition fTSl and . by (|2T)]) . Corollarv[T71 
could be extended to this kinds of random closed sets. 
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